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Abstract. When planning new facilities, one is interested 
to know whether and by how much the new technique 
is superior to already existing ones. We describe a gen- 
eral approach permitting to evaluate the relative merits 
of various techniques used in astrophysical observations, 
following the theory of model parameter estimation. It is 
applied to compare two-aperture optical long baseline in- 
terferometry to classical spectroscopy, both used to define 
the model parameters of the P Cyg wind. The wind mod- 
eling was done using an efficient approximation for com- 
putation of the line source function; it allowed to analyse 
about 10^ points in the parametric space of P Cyg en- 
velope models. It is shown that interferometry offers no 
real advantage if the object can be described by station- 
ary spherically symmetric models with a priori fixed ther- 
mal structure. However, if the object must be described 
by a model with a large number of free parameters, e.g. 
when the thermal structure of the envelope is not fixed a 
priori, then the interferometric measurements can reduce 
the error in model parameters determination by an order 
of magnitude. In the case of P Cyg, the reduction of er- 
rors provided by interferometry is highest for the baseline 
lengths in the range 45-90 m. This illustrates the capacity 
of the proposed method to be used for optimization of in- 
terferometric configurations. The simplifications adopted 
for this first trial application are critically analyzed, and 
future improvements are indicated. 
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1. Introduction 

When evaluating justification for new facilities, one may 
question whether and by how much the new technique 
is superior to already existing ones. We describe a gen- 



eral approach which leads to a quantitative answer, ex- 
pressed as a relative figure of merit of two observational 
techniques. The approach is based on the comparison of 
the corresponding errors in the model parameters deter- 
mination. The derived figure of merit can also be used in 
the closely related problem of optimization of complex ob- 
servations, thus contributing in maximizing their scientific 
output. 

We illustrate the developed method by comparing the 
technique of optical long baseline interferometry to the 
classical spectroscopy, both used to find the parameters 
of the gaseous envelope around the P Cyg star. For this 
first trial application, we use simulated data, obtained by 
solving the radiative transfer problem for grids of physical 
models of the envelope and then adding the photon noise. 

The chosen example of the interferometry vs spec- 
troscopy comparison appears to be timely. Indeed, in the 
recent years, interferometry has been undergoing a steady 
progress witnessed by publication of fringe visibilities for 
circumstellar envelopes of various types sometimes com- 
bined with a goo d spec t ral r esolution (Mourard etal 
19891 Vakil i et al. |l£ 

1994, |19971; Harmanec et al. |1996|; Stee et al 



1997|; Qu irrenbach et al 

19951 : 



1993 



also expected performances of interferometric arrays un- 
der construction, VLTI, von der Liihe et al. 1997, Petrov 
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et al. |1994 and CHARA, McAhster et al. |1997D . While 
this effort is fully justified, for the knowledge of the flux 
distribution on small angular scales can unveil important 
new features of astrophysical objects, it is also justified to 
wonder whether this effort is well payed for all classes of 
objects and for all kinds of physical problems. As a matter 
of fact, interferometric observations as compared to those 
at a single telescope are more complex, more expensive, 
more time consuming, and, when possible to compare, 
have a larger accumulated measurement error due to their 
complexity. Furthermore, in the foreseeable future the in- 
struments of this type will remain orders of magnitude 
less numerous and by all given reasons much less available 
for the community than single telescopes. Thus, it makes 
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sense to evaluate in quantitative terms what are the tar- 
gets for which the interferometric instrumentation has the 
highest potential in terms of astrophysical information as 
compared to single telescopes. 

Also, the optimization of interferometric observations 
may be a crucial aspect. Indeed, the existing and forth- 
coming long baseline interferometers comprise only few 
individual apertures, covering a small set of spatial fre- 
quencies during one observing cycle. On the other hand, 
the information provided by observations can depend crit- 
ically on spatial frequencies. Certainly, in such a situation 
the scientific yield of the interferometric array could be 
greatly improved if thorough model calculations can indi- 
cate the configurations expected to be the most useful for 
determination of aimed physical parameters of the object. 

Not surprisingly, the relative figure of merit derived in 
the present article depends not only on the experimental 
techniques, but also on the chosen set of model parame- 
ters, reflecting the fact that the physical model makes a 
necessary part of the merit problem. However, in modern 
astronomy, designing and building instrumentation on the 
one hand, and modeling objects on the other hand, were 
pushed to such a level of sophistication that they are as a 
rule two distinct activities. The optimum use of interfer- 
ometry needs a close collaboration of two communities. 

After these preliminary remarks, let us specify that 
in the present first trial application, our discussion will 
be limited to the comparison of observations with a two- 
aperture optical long baseline interferometer (hereafter 
OLBI) and the classical single telescope spectroscopy 
(hereafter shortly spectroscopy), both applied to study 
a circumstellar gaseous envelope formed by the wind of 
a massive star. The OLBI observations are assumed to 
have the same spectral resolution as spectroscopy. The 
discussion will be further limited to spectral profiles and 
visibilities in the Ha line of the P Cyg envelope, allowing 
to use published high quality observations both in OLBI 
and spectroscopic, and thus to provide a clear numerical 
illustration. 

In our earlier attempts to solve the OLBI vs spec- 



troscopy merits problem (Burgin & Chalabaev 1992; Bour- 
guine & Chalabaev 1994| ), we used qualitative comparison 
of observables, visibility and spectral profile, computed for 
a small number of envelope models. The results were am- 
biguous. Analising them, we arrived to the firm conviction 
that meaningful conclusions can be obtained only if (1) 
intercomparison of different types of observations is per- 
formed using a clearly defined figure of merit, indicating 
how much information on the studied object is provided 
by various techniques, and (2) a sufficiently large ranges 
in the space of envelope model parameters are analyzed. 

In the case of P Cyg, the required number of computed 
models (see Sect. ^) turned out to be 10** — 10^. An efficient 
simplified envelope model code, computing the emergent 
spectral profile and visibilities for the Ha line, was de- 
veloped and is described in Sect. 0. We considered only 



spherically symmetric outfiow models, the choice which is 
justified in Sect. ^. 

2. Relative figure of merit of observational 
methods 

2.1. Definitions and general theory. 

We consider astrophysical observations and their interpre- 
tation in the framework of the theory of parameter esti- 
mation (e.g. Kendall & Stuart 1967). That is, the observed 
object is assumed to be exactly described by a fixed mul- 
tiparametric physical model, the errors being solely those 
of the experimental measurements (see however Sect. 2.3). 
Any observable quantity is then a function of the vector 
of Mp model parameters = (Oi, 82, . . . , Oa/p)- 

An observation consists in finding the values of Mq ob- 
servables Y = {Yi , I2 , . . . , Ymo ) 1 where each observable Yi 
is either measured directly or can be calculated in a model- 
independent way as a known function of directly measur- 
able quantities. Obviously, observations with instruments 
of different types provide observables of different number 
and nature. 

The i-th observable can be represented as Yi — ^^(0) + 
£i, where Yi{&) is the "theoretical" value of correspond- 
ing observable in the absence of observational errors and 
e — (ei, £2, ■ • ■ J £j\/o) is the random vector of observational 
errors, which is assumed here to be distributed normally. 
We consider the case when Mq > Mp and the parameters 
could have been uniquely determined if the vector Y 
were known. 

We shall call the parameters to be determined the "tar- 
get parameters" and denote a set of Mt target parame- 
ters, or target set, as T. It is possible that Mt < Mp and 
the set T is a proper subset of the set of all parameters 
of the model. The situation of this kind may arise for two 
reasons. 

First, a model parameter may be set to some a priori 
fixed value because, say, it has been measured earlier by 
entirely different methods. This reduces the dimension of 
the problem and computational difficulties, the general 
method remaining the same. Let J- denote the set of Alp 
fixed parameters, and j be the value a priori assigned to 
a parameter Oy G JT. 

Secondly, some of the unknown physical parameters 
may eventually be considered "not interesting" if their 
values are not relevant to the astrophysical problem un- 
der study. Following the theory of parameter estimation, 
we will call them the "nuisance" parameters and denote 
the set of nuisance parameters, or "nuisance set" , as J\f, 
their number being Mn. As a rule, their infiuence can be 
separated in the error analysis only at the final stage of 
computations, for their values are to be calculated along 
with the values of target parameters. 

Interpretation of an observation consists in calculating 
the vector &t{Y) = (Ot(i), Ot(2), ■ ■ • , ©((Mt)) that pro- 
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vides the best fit to the observed values Y. Here t{i) is 
the index of the component of vector Q corresponding to 
the same model parameter as the i-th target parameter 
(z = 1,...,Mt). 

The precision of the values 0t resulting from interpre- 
tation of the observations is characterized by M (0t — ©t) 
and D0T, where for any random vector ^ expressions M£ 
and denote its mean and its covariance matrix respec- 



tively, and ©T = (0t(i), 



t(l),Ct(2),. 



t(MT 



In what follows, we assume that Y is an unbiased esti- 
mate of y , i.e. Me = 0, and that the measurement errors 
are small enough, so that the error analysis can be per- 
formed using the linearized version of the least squares 
method. 

We additionally assume that all a priori fixed param- 
eters are set to their true values (see however 2^), that is 
Qf = Qf for any 0/ G JF, then M0 = and statistical 
properties of errors in parameter determination are com- 
pletely characterized by the covariance matrix of errors C 
= D0. 

2. 2. Random errors and the relative figure of merit 



According to Kendall & Stuart ( 1967 , Chap. 19), the co- 
variance matrix of errors C is related to Y and by the 
following expression: 



c(0,r,AA) = (aT( 



(1) 



where A is the A/q x (A/t + A/n) matrix with elements 
Aij = for Qj e T U A/" and e is the vector of experi- 



mental errors defined in Sect. 2.1 



The natural scalar characteristics of the precision of 
parameter determination for a given set T of target pa- 
rameters is the following principal subdeterminant of the 
covariance matrix: 



C(0,T,A^) 



(2) 



det 



Q(i)t(i) Q(i)t(2) 

Ct(2)t(l) Q(2)t(2) 



-t(MT)t(l) Ct()t(2) 



-t(l)t(MT) 

-t(2)t(MT) 



-t(MT)t{MT) 



Geometrically, this subdeterminant is proportional to 
the hypervolumc of the scattering ellipsoid in the space 
of target parameters. It depends not only on the physical 
model and errors of measurements, but also on the ana- 
lytical form used for description of the model: two phys- 
ically equivalent but mathematically different (e.g. inter- 
related by a reversible substitution of variables) analyti- 
cal representations could yield entirely different values of 
C(0,r,AA). 

However, a pair of observational methods can well be 
compared if one makes use of the ratio 



R{@,T,M) = (Ci(0,r,A/')/C"(0,r,A/')) 



l/M-T 



where and C^^ are the subdeterminants calculated for 
the observation of the same object with instruments I 
and II. This ratio depends only on the physical model 
used, and the obtained precision on the model parame- 
ters. We shall call the quantity R{@,T,Af), introduced 
by the Eq. (^ , the "random error ratio" , or the "relative 
figure of merit" . 

If the instruments are of the same kind, and differ 
from each other only in precision, the value R{@,T,Af) 
is merely the ratio of observational errors. If the instru- 
ments are of different kinds, providing observables of dif- 
ferent nature and number, and the target set consists of 
only one model parameter, the value R{&, T, M) is merely 
the ratio of resulting random errors in the parameter de- 
termination. However, in the general case of instruments 
of arbitrary kinds and multiparametric models, no simple 
ratio of errors exists, and the evaluation of relative merits 
can be done only using the quantity R{&,T,J\f) defined 
by Eq. (|). 

It allows to compare various observational techniques 
applied to objects described by various multiparametric 
models: equality R{&,T ,Af) < 1 means that the instru- 
ment I is better suited for determination of parameters 
from the target set than the instrument II. 

2.3. Systematic errors induced by interpretation 

The fact that the description of an object by a multipara- 
metric physical model is only an approximation to the 
reality implies that we have to consider the robustness of 
the method, that is the stability of the results it yields 
with respect to deviations of the real situation from the 
model. 

The present framework offers a way to obtain cer- 
tain quantitative characteristics of the robustness. Indeed, 
let us consider a multiparametric model with ^ 0. If 
9/ 7^ 6/ for some Qf G T then, in general, MG)t — 9f ^0 
for Qt G T. That is, in addition to random errors of ob- 
servational origin, the result is biased by systematic errors 
due to inaccurate interpretation. The value of that bias 
characterizes the robustness of the method with respect 
to deviations of 9/ from its true value. 

In the linear approximation, M(9t — 9() = S x (9/ — 
9/), where S is the Mt x A/p matrix with elements 



99 



/ 



(3) 



where 9t e r, 9/ G J^. 

A comprehensive study of systematic errors requires a 
joint analysis of individual elements of the matrix S and 
uncertainties in model parameters from J-. In the present 
paper, we will develop a simplified approach providing 
semiquantitative indications concerning the relative ro- 
bustness of different observational techniques. 
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Let us first define for each observational method under 
consideration the value 



c/(e,T,^) = (c(e,r,^)/c(e,r,( 



nI/A/t 



(4) 



When the estimates for target parameters are uncorrelated 
with estimates for parameters from T ^ equality U — \ 
takes place: the method is robust (with respect to devia- 
tions of specified form from the model). When such a cor- 
relation exists, V exceeds unity; the closer the correlation, 
the larger its value, eventually implying poor robustness 
of the method. 

Further, to compare the robustness of methods I and 
II with respect to inaccuracies in parameters from JF, we 
will define the "robustness ratio S"' as follows: 

5'(0,T,j^) = J/VC/" = i?(0,r,J^)/i?(0,r,0) (5) 

where function R is defined in Eq. (^. Note that once 
the values of -R(0, T, A/"), which are necessary for analysis 
of random errors, are found, comparison of robustness is 
straightforward. 

3. The relative figure of merit of OLBI vs 
spectroscopy 

In this section we obtain the expressions for the elements 
of the covariance matrix D© in the form that can be di- 
rectly used to compare the OLBI and spectroscopy. 

First, the measurement errors have to be specified. For 
this first trial application, the only considered source of 
random errors is the photon noise. Other sources of er- 
ror, in particular those in OLBI arising from atmospheric 
seeing and calibrations of the modulation transfer func- 
tion (see Roddier & Lena 1984 ; Mourard et al. 1994 ) are 
important and has to be incorporated in a future work. 

Also, the spectral coverage and the spectral resolution 
have to be specified. We assume that they are the same 
for both kinds of observations; the covered spectral region 
is [Ai,A2], the spectral resolution is sufficiently high, i.e. 
the monochromatic intensity received from each point of 
the observed object remains nearly constant across each 
spectral channel. In other words, we assume that AA, the 
width of the spectral channels, in terms of velocity is lesser 
than the thermal velocity of the emitting atoms. For the 
P Cyg type envelopes this corresponds to AA ^ 0.3 A, 
which is easily satisfied with modern spectrographs. Then 
one can formally consider the convenient limit case AA — > 
0. In this case, the resulting covariance matrices tend to a 
limit which is independent of the spectral resolution (cf. 
Eqs. (|ll|) and (|l|)). 

3.1. Spectroscopy 

The physical quantity provided by spectroscopy is the flux 
density as a function of wavelength. In terms of observ- 
ables, it is given by the vector with the components 

f,s = ^„ i = l,2,...,L, (6) 



where Ni is the number of photons recorded in the i-th 
spectral channel during the exposure, and L is the num- 
ber of channels. The random values Ni are related to the 
physical parameters of the object by 



where 

N,^F{@,\,)EA,X, 



(7) 



(8) 



F(8, A) is the monochromatic flux from the object at 
wavelength A, A^A is the width of i-th spectral channel, 
and Tii is the measurement error due to the photon noise, 
with the variance given by 



Dni ^ Ni . 

Finally the factor E is given by 
E = Acst , 



(9) 



where Acs is the instrument effective aperture and t is the 
exposure time. Hereafter we suppose that = E^, where 
E^ and E^ pertain correspondingly to the spectrometer 
and interferometer. Since the values E^ and E^ enter the 
final results only through the ratio E^ /E^, we can for the 
sake of simplicity set E^ — E^ = 1. 



-S _ 



trix of measurement errors are given by 



1(ef4) 



ifi^j 
iV, if ^ = j 



(10) 



Substituting Eqs. (|6[)-(p^) into Eq. (|^) and proceeding to 
the limit A^A ^ we obtain: 



' pq 



1 dF{&,X)dF{&,X) 



^ F(0,A) dOp 



dX. 



(11) 



For spherically symmetric objects considered in the 
present paper. 



/•OO 

F{e,X)=27T I{Q,X,p)dp, 
Jo 



(12) 



where /(0,A,p) is the monochromatic intensity at the 
angular distance p from the center of the object. 

3.2. OLBI 

We consider the case of a two-aperture OLBI with an ad- 
justable baseline. The physical quantity eventually pro- 
vided by OLBI measurements is the fringe visibility V 
as a function of the baseline vector. The visibility V is 
equal to the real part of the complex degree of coherence, 
and therefore is the Fourier transform of the brightness 
distribution in the focal plane by virtue of the Van Cit- 
tert - Zernic ke theorem (e. g. Mar iotti 1998, textbooks of 



Pefina |1972| , and Goodman [19851) . However, the value es- 
timated in practice is V^. We refer the reader for details 
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to the thorough discussion by Mourard et al. (1994). In 
the present analysis we consider the estimates Wi of the 
quantity 



W, = W{@, \,B) = V^{@, X„ B), 



(13) 



where i refers as previously to the i-th spectral channel, 
and B is the projected baseline length of the interferom- 
eter. In the case of spherically symmetric objects the in- 
tensity distribution is circularly symmetric and V is given 
by th e normalized Hankel transform as follows (Bracewell 



1978 



V{&,X,B) = 



2tt I{@,X,p)Mkp)dp 

F(e,A) 



(14) 



where Jo is the Bessel function, and k = 2'kB/\. 

Statistical errors Wi = Wi — Wi obey the equations 



M{wi'Wj) = 
and 

U{winj) = 







iWi (r, 



if « ^ j 



(15) 



(16) 



The validity of Eqs. ( |15|) a nd ( |16| ) for i ^ j is evident. The 
case i = J is treated in lAl. 

As a rule, the spectrum of the object, i.e. the values Ni, 
is also recorded. Consequently, the vector of observables 
for the OLBI and the elements of its covariance matrix of 
errors are given by 



and 





■,Nl,Wl 



(17) 



if « ^ j 
if i = j = 2k 
^{2-Wk) iii = j = 2k 



Substituting this expressions into Eq. (^, we obtain in the 
limit AiA ^ that 



Ipq 



(C«) 



PI 



F 



dW dW 



i^^ AW{2-W)dQpdQq 
where (O?)^^ is defined in Eq. (|l 



dA, (18) 



3. 3. The relative figure of merit 

As discussed in Sect. |2.2| , the relative figure of merit of 
OLBI as compared to spectroscopy is the random error 
ratio: 

R^^{&,T,Af,B) ^C\&,T,Af,B)/C^{&,T,M). (19) 

From Eqs. (|l|) and (|l|), it follows that R^^ < 1, which 
merely reflects the fact that the observables of spec- 
troscopy constitute a subset of the OLBI observables. 
Consequently, in the adopted comparison, the statistical 



errors in parameter determination by OLBI are always 
lesser than those obtained when only spectroscopic data 
are used, if however the spectroscopic data are of equal 
precision. 



As explained in Sect. 3T, the Eq. (19) is obtained in 
the assumption that = . To compare a pair of in- 
struments for which that equality does not hold, the RHS 
of Eq. dH) should be multiplied by {E^/E^f/^. 



4. Models of the P Cyg v^rind. Computed 
observables 

4-1. Stellar parameters 

The star P Cyg (Bl lape) is characterized by a high mass- 
loss rate and exhibits bright optical emission lines, formed 
in the dense and nearly fully ionized stellar wind. It be- 
longs to the class of Luminous Blue Variables, a short- 
living transition phase in the evolution of a massive star 
at the end of the hydrogen burning (see e.g. Humphreys 
& Davidson |1994 Langer et al. |1994| , Maeder |1997| ). 

There were at least two strong reasons to choose P Cyg 
as the astrophysical object of the present study. First, 
this bright object, V = 4.8™, is one of the best stud- 
ied emission line stars. Although in this first trial we use 
simulated data, it is important to note that for P Cyg 
there exists a rich literature providing not only high qual- 
ity spectroscopic data (e.g. Scude ri et a l. 1994 ), but also 
interferometric data (Vakili et al. 1997), as well as thor- 
ough theoretical analysis (Drew 1985; Pauldrach & Puis 
119901 ). This insures that the present work can be followed 
up by a practical application. Secondly, it happens that 
in the growing list of outflows known to be non-spherical 
(e.g. Wolf, Stahl, FuUerton |1998D , that of P Cyg is an ex- 
ception exhibiting spherical symmetry to a good degree 
of accuracy (Nota 1998), becoming clumpy only on short 
time and small flux scales (Taylor et al. |1991l Vakili et 



al. 1997; Nota 1998] ). This implies that the assumption of 
spherical symmetry of the envelope, used in the present 
work, is realistic. Let us remind that it allows to com- 
pute a model in a reasonable amount of time, and thus 
to explore a large domain of parametric space, which is a 
requirement for the evaluation results to be meaningful. 

For the distance to the star d, its radius i?* and the 
effective temperature T^a we adopted the following values: 
d = 1 800 pc, i?, = 76.0 j ^0, Tc ff = 20000 K (Lamers et al. 
1983| Pauldrach & Puis |1990| ). The spectrum of the star 
was assumed to be blackbody. Our study is limited to the 
hydrogen Ha line, which is the most prominent and the 
best studied feature in the spectrum of the star. 

4. 2. The physical model of the envelope 

The mass distribution within a spherically symmetric and 
stationary outflow is described by its mass loss rate M 
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and the outflow velocity v, which we assume to depend on 
radial distance r in the following way: 

v{r)=Vc + {voo-Vc){l^l/r)", 

where Vao is the terminal velocity of the wind, is ve- 
locity of the wind at the base of the photosphere, and a 
is a dimensionless parameter that characterizes the rate 
at which the velocity approachs its asymptotic value at 
large r. The radial distance r is measured in units of stel- 
lar radius i?,, while M, Vc, Woo, and a are parameters of 
a model. 

The temperature in the envelope is often assum ed to be 
constant. However, as it was shown by Drew ( 1985 ), across 
the region of the Ha line formation (i? ^ 5) the tempera- 
ture in the envelope can vary by as much as 6000 K. There- 
fore, along with isothermal models, we also computed the 
emergent emission for non-isothermal models. 

To keep a finite number of scalar parameters, the de- 
pendence of envelope temperature T{r) on radial distance 
is approximated in the following manner: it is assumed 
that T{r,) = Tj for ri = 1 < r2 < . . . < rjif^, that T(r) 
is a linear function of logr at each interval [ri,ri4-i] for 
< i < Mi{ - 1, and that T{r) = Tm^ for r > ta/^. The 
values Ti were fixed for each model of the family. 

We consider only non-increasing temperature laws 
T(r), that is A,T > for i = 2, . . . , Mr, where A^T = 
Ti^i — Ti. Since the computations are organized in such 
a way that the model parameters vary independently on 
each other, it is the values A^T that are used along with 
Ti as the parameters defining the envelope temperature, 
and the vector of model parameters is given by 

= (a, M, Voo,Vc, Ti, A2r, . . . , AmrT). 

4-.3. Source function and Radiative transfer equation 

The "supersonic" Sobolev approximation is a well known 
efficient method of calculating line spectral profiles (e.g. 
Castor f970 ). However, it would have a low accuracy if 
applied in the P Cyg case, for a noticeable fraction of 
Ha flux is emitted in the inner part of envelope, where 
the outflow velocity is comparable to that of sound. We 
adopted a mixed method computing the source function 
in the Sobolev approximation, and then finding the emer- 
gent intensities by exact numerical solution of the transfer 



equation. As it was shown by Hamann (f981), the errors 
in line profiles computed in the original Sobolev approx- 
imation come mainly from calculations of the emergent 
intensities, whereas the source function is accurate for a 
wide range of physical conditions. 

The choice of the model of the hydrogen atom was 
based on the fact that the regions of the envelope emitting 
the major fraction of the Ha flux are nearly completely 
ionized and opaque to Lyman continuum and La, so that 
direct recombinations to and photoionizations fro m th e 
ground level cancel out, and La is saturated (Drew 1985|) . 



The populations of levels n = 2 and n = 3 are mainly 
defined by coUisional and radiative transitions between 
these two levels, radiative ionizations due to stellar radia- 
tion, and radiative recombinations (including indirect) to 
the level n = 3. We therefore adopted the three- level + 
continuum model of hydrogen atom. 

The balance equations for 712 and 713 , the number den- 
sities of hydrogen atoms respectively at levels 2 and 3, 
take then the following form: 



"•2 (C23 
-"2 C23 



+ I2) 



/332^32) — R2 

+ I3) — R3 



(20) 



"-3 ((^32 
"-3 (C32 

where ^32 is the spontaneous emission coefficient, C23 and 
C32 are collisional excitation and de-excitation coefficients 
for the transition, /2,/3 and i?2,^3 are respectively ion- 
izations coefficients and recombination rates for the levels 
involved, and (332 is the escape probability. 

Since the hydrogen is nearly completely ionized, the 
values C23, C32, R2 and R3 can be considered independent 



on n2 and n^, the only source of non-linearity in Eqs. (20) 
being the terms containing /323- 

The escape probability P32 is a rather complex func- 
tion of the population of the lower level of the transition, 
so that commonly the Eqs. (20) are solved by iterations, 
and this and, in particular, the recalculation of the escape 
probability, is by far the most time consuming operation. 

To speed up computations, we developed a new ap- 
proximate method for solving the equations of statistical 
equilibrium. It is based on the fact that, as shown in ^ 
the function ^32(^2) can be approximated with sufficiently 
good accuracy by a simple analytic expression as follows: 



/332(»^2) = 1/(1 + n2/n^s), 
where 

Stt 92 / 1 dw(r) 2v{r) 



(21) 



(22) 



>?M2 33 \3 dr 3 r 

(72 and 53 are statistical weights of the levels, and A is the 
line wavelength. 

Introducing the dimensionless variables 



P2 = n2/"-o : 
where 

i?2 

"0 = -i-, 



P3 = n-s/inoQ) , 



(23) 



23 



93 

— cxp . 

92 \ kT 



E. 



and £'23 is the energy of transition, and substituting the 
approximation (21) into Eqs. (pO[) we obtain after some al- 



gebra that P2 is the positive solution of the quadric equa- 
tion 

aP| -(- &P2 + c = , 
where 

a = -~p2MS, 

b = (p2*(l-/3o*('5 + 5)) + Af *(f-/3o)+P3)^, 
c = /3o * {iP2 + P3) *{l + S) + M) = . 
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Here 

Pi = ^i/(^32»T-o) for i = 2,3, 
M = lpIS, 
S = 1 + lQ, 

f3o = 1/(1 + no/^as) . 

and 

i = h/h ■ 

When P2 is calculated, the value of P3 can be found using 
the relation 



ft = 



QP2 + p:i5 
Q(lp25 + P5+1) 



The source function is then easily computed from level 
populations, which can be obtained using Eq. (|2^). 

Integration of the transfer equation was performed us- 
ing the code developed by Bertout(1984), who kindly pro- 
vided it to the authors. In computing the line profile, the 
code assumes the envelope to be isothermal, so that its 
use for a non-isothermal case requires some comments. 
The envelope temperature enters the calculations at two 
points: (1) In calculating the source function, through co- 
efficients C23 and C32 of Eq. (20). Since the Bertout code 
is applicable to arbitrary source function, variations in 
temperature does not cause any difficulties here. (2) In 
integrating the transfer equation, T(r) enters the result 
through the local Doppler line width. Since we consider 
only the envelopes with relatively low temperature con- 
trast (Ti — Tmr)/Ii ^ 0.4, and dependence of the Doppler 
width on temperature is rather week, the error induced by 
this isothermal code structure is still negligible. 

4.4. Computation of final results 

Once the intensity /(0, A,p) is computed, the observable 
quantities, i.e. the emergent flux i^(0, A) and the visibility 
squared W{@,X,B) are obtained using Eqs. (|l^)-(|l^. 

The relative figure of merit i?(0m,T,A/',-Bj) and the 
robustness ratio S{&„i,T,J^,Bj) are calculated in two 
stages, implemented as separate programs. 

First, for the given grids of model parameters {©,„ : 
m — 1, . . . , Nyn} and the projected baseline lengths {Bj : 
j — 1 .... , Nb } , we calculate and store the matrices 



^ and (C^) \ defined in Eqs. 1^ and (|T|). The 



defining parameters of the grids are entered as input data, 
the grid of models being constructed as the direct product 
of uniform grids for individual parameters. 

The derivatives entering definitions of the matrices are 
approximated by finite differences. All the model param- 
eters that are not constant on the grid are treated as ad- 
justable, that is either target or nuisance. The distinction 
between those two types of parameters is irrelevant at this 
stage. 



At the second stage, for given partitions of set 
T U A/" of all adjustable parameters on the subsets T 
and M, we compute the values of R{&„i,T,JV,Bj) and 
S{@„^,T,T,Bj) (see Eqs. (|) and (|)) on the grid. In this 
way, the dependence of the results on T and JV can be 
studied without repeating the time consuming physical 
modeling of the envelope. 

5. Results 

In total, the values R{Q, T,Af, B) and S{@, T, B) were 
calculated as a function of projected baseline length for 
various sets of target and nuisance parameters at more 
than 10^ points of the parametric space. Table |l] presents 
the principal characteristics of the model grids. The ranges 
of variations of model parameters on the grids were chosen 
so as to include the values of parameters from earlier works 



on P Cyg cited in Sect. 4.1 



The main result of the present study is that the in- 
crease in accuracy provided by OLBl strongly depends on 
the number of free parameters in the model: it varies from 
almost negligible for models of low dimension to very sig- 
nificant for models of high dimension. Since in our case the 
boundary between "low" and "high" is nearly coincident 
with division onto models with fixed and adjustable ther- 
mal structure, we discuss the results pertaining to those 
two classes of models separately. 

5.1. Models with fixed thermal structure 

Among our grids, the grid 1 most closely corresponds to 
the models often used in interpreting Ha observations of 
P Cyg stars, where foo is excluded from the set of ad- 
justable parameters of the model (see e.g. Scuderi et al. 



1994) and its value is taken from the analysis of absorp- 



tion lines of metals in the ultraviolet (Casatella et al. 1979 
Lamers et al 1985). 



The typical results for an individual point in the para- 
metric space are illustrated by Fig. where the ran- 
dom error ratio R^^ = C^{@,T ,M . B) /C^ {@,T ,M , B) 
is plotted as a function of the normalized baseline length 
B/Bq for various T and Af. Here B is the projected base- 
line length, Bq — \/2Tr6, and S is the angular radius of 
the central star. For the adopted values of i?* — 76 Rq, 
d = 1800 pc, and A = 6265 A one gets Bo = HOm. 

The principal features of dependence of R^^ on B are 
independent on details of physical model used and can 
be easily explained qualitatively. As shown in Sect. |3.3| , 
inequality R^^ < 1 holds for all < B < 00. When 
the projected baseline length is very small, the interfer- 
ometry does not provide any additional information as 
compared with the spectroscopy, for the object gets unre- 
solved. Consequently, R^^ ^ 1 for B ^ 0. In the opposite 
case of large B the fringe contrast is close to for all phys- 
ically realistic intensity distributions and its dependence 
on model parameters is impossible to measure because of 
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Table 1. Grids of models: the number of points in the parametric space N^; dimension, that is the number of varied 
parameters, of the grid D] the values of parameters that were fixed on the grid; and the ranges of variation for variable 
parameters. 



Grid 
no. 




D 


a 


M 


km s-^ 


Vc 

km s-^ 


Ti 
K 


i?2 


K 




K 


i?4 


A4T 
K 




AsT 
K 


1 


10 000 


3 


3.5 
4.5 


1.010-5 
2.010-5 


200 


15.0 
55.0 


11000 
15 000 


















2 


10 000 


3 


3.5 
4.5 


1.010-5 
2.010-5 


200 


15.0 
55.0 


11000 
15 000 


2.0 


1300 


3.4 


2100 


5.0 


200 


10.0 


3000 


3 


32 805 


4 


3.5 
4.5 


1.010^5 

2.010-5 


150 
250 


15.0 
55.0 


11000 
15 000 


















4 


15 625 


6 


3.5 
4.5 


1.010-5 
2.010-5 


200 


15.0 
55.0 


11000 
15 000 


2.0 




1300 


5.0 




2300 










5 


46 656 


6 


3.25 
4.25 


1.010-5 
2.010-5 


200 


15.0 
55.0 


11000 
15 000 


2.0 




1300 


3.4 




2100 










6 


6561 


8 


3.5 
4.5 


1.010-5 
2.010-5 


200 


15.0 
55.0 


11000 
15 000 


2.0 




1300 


3.4 




2100 


5.0 




200 


10.0 




3000 


7 


65536 


8 


3.5 
4.5 


1.010-5 
2.010-5 


200 


15.0 
55.0 


11000 
15 000 


2.0 




1300 


3.0 




2100 


4.0 




200 


5.0 




3000 


8 


6561 


8 


4.0 
4.5 


1.510-5 
2.310-5 


180 
220 


20.0 


14000 
16 000 


2.0 




1000 


3.0 




1000 


4.0 




1000 


5.0 




1000 


9 


19 683 


9 


3.5 
4.5 


1.510-5 
2.310-5 


220 
300 


15.0 
25.0 


16 000 
18 000 


2.0 




1300 


3.0 




2100 


4.0 




200 


5.0 




3000 




0.8 ' ' ' ' 1 

O 0.5 1 1.5 2 

Normalized base length 

Fig. 1. Error ratio for isothermal model with a = 3.94, 
M = 1.441O-^M0yr-i, = 32.8 km s'^, Ti = 12778 K, 
Voo — 200 km s-^, and the set of adjustable parameters 
TuAf = {a,M,Vc,Ti}. Target set T = {a} (solid line), 
{M} (dashed hue), {Ti} (short-dashed line), and {vc} 
(dotted line). 



observational errors. That is in that case the interferome- 
try again provides no additional information and R^^ 1 
for B ^ 00. 

It follows from the above stated general properties of 
R^^ that for any T and A/" it has a (possibly non-unique) 
minimum at a certain finite value of B. The location of 
the minimum indicates the projected baseline length at 
which interferometric observations are most informative, 
and the corresponding value of R^^ characterize to what 



degree the interferometry can reduce the errors in model 
parameters determination as compared to spectroscopy. 

Let us note that it would be insecure to compare the 
merits of the methods using values of R^^ obtained for 
an isolated point or even for few arbitrarily chosen points 
of parametric space. Since the model parameters are not 
known a priori, it is necessary to use integral characteris- 
tics describing the behavior of R^^ on the whole grid. We 
will use the following three values calculated as a function 
of Bj/Bq for various T and Af: 

i?„,i„(T,A/-,S,) = min R'^i&„^,T,M, B,) (24) 

1 < m < Nm 

i?max(T,AA,S,) = max R''^{&„„T,M,B,) (25) 

1 < m < Nj^ 

Popt(r, A/-, B,) = A^opt/A^m, (26) 

where A'opt is the number of points on the grid for which 
function R^^{Q, T,J\f, B) reaches its minimum bX B = Bj. 

The function defined in Eq. ( p6| ) is immediately related 
to the problem of optimal choice of the baseline length: for 
given T and A/", the higher the value Popt{T,M, B), the 
higher the probability that the interferometric observa- 
tions at the baseline length B would yield most accurate 
model parameters. 

The dependence of i?min and Popt on the normalized 
baseline length and set of target parameters for grid 1 is 
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in 
d 



in 
d 



d 



in 
d 




/ 



H 1 1 1 I 1 1 h 



H 1 1 1- 



iZ=i L 



0.5 



1.5 



Normalized base length 

Fig. 2. The lower limit of random error ratio -Rmin (solid 
line) and distribution of optimal baseline lengths Popt 
(dashed lines, plotted in an arbitrary scale) for grid 1. 
The set of adjustable parameters T U A/" = {a, M,Vc, Ti}, 
target parameters are indicated on the plots. 




E 

d 



d 



/ 



\ 



H h 



I I I I I I 



-I \ h 



a,Vc,Ti, A2T, A^T 



0.5 



1.5 



Normalized base length 

Fig. 3. The lower and upper limits of random error ratio 
(solid lines, i?max differs noticeably from unity only in the 
bottom plot) and distribution of optimal baseline lengths 
Popt (dashed lines, plotted in an arbitrary scale) for grid 5. 
Target set T = {M}, nuisance parameters are indicated 
on the plots. 



shown in Fig. |2[ For this grid, the value of i?max is close 
to unity for all T and Af and is not shown in the plot. It 
can be easily seen that for this grid of models OLBI pro- 
vides only relatively small reduction of the random error in 
parameter determination as compared with spectroscopy. 
This conclusion remains also valid if we add the terminal 
velocity Voo to the set of adjustable parameters (grid 3). 

As i?^^(0, T,J\f, B) « 1 for aU T and TV, the same ap- 
proximate equality is evidently valid for the robustness ra- 
tio S{@,T,T,B) = R{@,T,T,B)/R{@,T,i/},B), which 
indicate (see Sect. ^.3[ ) that, for models with a priori fixed 
thermal structure, interferometric data will neither reduce 
significantly the level of systematic errors. 

5.2. Models with adjustable thermal structure 

The number of possible target sets is an exponentially 
increasing function of the dimension D of the model pa- 
rameter space, so that when the values A^T are added to 
the set of unknowns, it becomes impossible to make an ex- 
haustive presentation of even the integral results, i.e. the 

functions i?rnin(^, A/", B), 

R 

max 

{T,Af,B), and Popt (T,AA,P). 
In Fig.||, we show the dependence of Pmin, Pmax, and 
Popt on the number of nuisance parameters for one of our 
simplest grids of models with adjustable thermal struc- 
ture, the target set consisting of one parameter, the mass 
loss rate M. From many points of view, it is the most im- 
portant parameter of the stellar wind, characterizing both 
its effect on the evolution of the star and the influence 



of the outflowing matter on the surrounding interstellar 
medium. 

As it can be seen, the lower limit of the random error 
ratio Pjnin rapidly decreases as the number of adjustable 
temperature parameters increases, so that when the ther- 
mal structure of the envelope becomes a nuisance param- 
eter, interferometry can reduce the random error in M by 
as much as an order of magnitude. 

Note, that the decrease in Pmin is caused by the fact 
that parameters A^T are adjustable, rather then by devi- 
ation from isothermicity. This can be seen from the upper 
two plots in Fig. where Pmin > 0.8 for a grid contain- 
ing high proportion of models with noticeable temperature 
gradient, and is also supported by the results obtained for 
grid 2, which differs from grid 1 only in that its thermal 
structure approximates that of Drew's model B, which 
corresponds to a = 4, t;c = 15km/s, v^x, — 300km/s, and 
M = 1.510^ A/Qyr^^ in our notations, instead of being 
isothermal. 

For other single-parameter target sets the dependen- 
cies of Pmin and Pmax On M are qualitatively the same. 
Fig. ^ displays the results for the most important case 
when all the parameters of the model except for v^o are ad- 
justable. As one can see, again, in favorable circumstances 
interferometry can increase accuracy by an order of mag- 
nitude. In contrast to Pmin, the value i?max is close to 
unity for all projected baseline lengths and combinations 
of sets T and J\f that we have studied. This signifies that 
for any projected baselength there exists a combination of 
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Normalized base length 

Fig. 4. The lower and upper limits of random error ratio 
(solid lines) and distribution of optimal baseline lengths 
Popt (dashed lines) for grid 5, T U7V={a, M, Wc, Ti, AaT, 
A3T}, and target sets indicated on the plots. 



model parameters for which interferometry practically do 
not provide reduction in error. 

It is interesting to notice that for target sets containing 
several parameters (see Fig. ^) , the influence of individual 
parameters on i?inin and i?max to some extent averages 
out: i?min IS Systematically higher, and i?max is systemat- 
ically lower than the corresponding values for individual 
parameters. For this reason, when more than one model 
parameter is to be determined, interferometry at a single 
projected baselength is unlikely to yield overall gain in ac- 
curacy in excess of factor of two even at optimal baselines. 

We studied further the effect of increasing the number 
of sublayers Mr, which permits to investigate the influ- 
ence of fine details of the thermal structure of the enve- 



lope. This appears to be justified for Drew (1985) showed 
that in certain cases this structure is quite complex, and 
gradient of log T(r) varies strongly with radial distance. 
Consequently, a model aiming to closely approximate the 
realistic envelope should be parameterized using a large 
Mr (see Sect. 4.1) and hence requires computations on 



grids of high dimension D. 

Since the number of mesh points of the grid depends on 
D exponentially, for computations based on such complex 
models (grids 5-9 from Table |I]) the limited computer re- 
sources dictate sparser mesh points along each axis of the 
parameter space and a tuning of parameters controlling 
the computational process in such a way as to increase 
the speed of computation at the expense of precision. 




J 1 1 1 1 I 1 1 1 i_ 



Normalized base length 

Fig. 5. The lower and upper limits of random error ratio 
(solid lines) and distribution of optimal baseline lengths 
Popt (dashed hues) for grid 5, T UA/'={a, M, v^, Ti, A2T, 
AsT}, and target sets comprising several parameters (in- 
dicated on the plots). 

Although numerically less precise than computations 
for simpler models, our results for more realistic grids indi- 
cate that the conclusions derived for low Mr can be safely 
extrapolated for higher Mr. This point is illustrated in 
Fig. ^ showing the results for single-parameter target sets 
for grid 7. As it can be seen in Figs. |3|-|| and is confirmed 
by our other results not presented here, the distribution 
of optimal projected baseline lengths in the majority of 
cases has a maximum in the range 0.4-0.8 in dimensionless 
units. In some cases this maximum shifts to higher base- 
line lengths, but this occurs only for a very special choice 
of T that correspond to measurements of wind tempera- 
ture gradient at the distances of several stellar radii from 
the star with all other model parameters considered nui- 
sance (see e.g. bottom plots in Figs. ^ and Thus, the 
aforementioned range, corresponding to 45-90 m in linear 
scale, can be recommended as the optimal choice for inter- 
ferometric observations aimed at determination of global 
characteristics of the P Cyg wind. 

The typical dependence of the robustness ratio 
S^^{@,T,M,B) = U^/U^ on the baseline length and the 
set of fixed parameters for a model with adjustable ther- 
mal structure is shown in Fig. M. At nearly all points of 
parametric space S^^ < 1, and iS is generally lower when 
the set of fixed parameters !F includes the values A^T. For 
a small fraction of models and for short projected baselines 
where interferometry can not significantly reduce random 
error, the value S^^ > 1 by a negligible amount. 
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Normalized base length 

Fig. 6. The lower and upper limits of random error ratio 
(solid lines) and distribution of optimal baseline lengths 
Popt (dashed lines) for grid 7, T UAf={a, M, Vc, Ti, A2T, 
A3T, /S.4T, A5T}, and target sets consisting of a single 
parameter (indicated on the plots). 




Normalized base length 

Fig. 7. Robustness ratio for the model a = 4.00, M = 
1.5 10-5 Mg/yr, Vc = 35.0 km s'^, Ti = 13000K, A2T = 

650 K, A3T = 1050 K, AaT = 100 K, A5T = 1500 K from 
grid 6 and T={M}. Set of fixed parameters T = {a, Wc, 
Ti, A2T} (solid line), {a, Vc, Ti, A2T, A3T} (dashed line), 
{a, Uc, Ti, A2T, A3T, A4T} (dash-dotted line), and {a, 
Vc,T-L, A2T, A3T, A4T, A5T} (dotted line). 

Thus, our results indicate that if the baseline length 
and the model used for interpretation are chosen in such 
a way as to reduce random errors, OLBI appears also able 
to reduce the systematic errors. 
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6. Conclusions 

Using the theory of model parameter determination, we 
developed a general method which permits a quantita- 
tive comparison of observational techniques and optimiza- 
tion of complex observations through the relative figure of 
merit, defined as a generalization of the ratio of random 
errors of model parameters. 

The method was applied to compare the outcome of 
OLBI and classical spectroscopic observations, both used 
to determine the parameters of the outflow of the P Cyg 
star. The observable quantities were computed using an 
efficient radiative transfer code and realistic grids of the 
envelope model. The OLBI and spectroscopic measure- 
ments errors were simulated assuming that the only source 
of errors is the counting statistics, i.e. we considered a 
nearly-ideal interferometer. Other main assumptions and 
simplifications of the present work were as follows: (1) the 
P Cyg wind is stationary and spherically symmetric; (2) 
the only OLBI observable considered here is the visibility 
modulus; (3) only the hydrogen Ha line is considered. 

The main conclusions of the present work are as fol- 
lows: 

1. A meaningful evaluation of the OLBI vs spectroscopy 
relative figure of merit requires exploration of large do- 
main of the model parametric space (about 10^ points 
computed in the present work). 

2. If the P Cyg outflow can be accurately described by 
a simple model of a stationary spherically symmet- 
ric isothermal envelope, then the use of interferometry 
does not substantially reduce the random error in the 
model parameter determination. 

3. If the P Cyg outflow requires more complex models, 
e.g. with a priori unknown dependence of the envelope 
temperature on radial distance, use of a nearly-ideal 
interferometer can reduce errors in model parameters 
by an order of magnitude as compared to the spec- 
trometer of the same collecting power and noise level. 

4. The study of the relative robustness, i.e. the stability 
of the obtained model parameter determination with 
respect to systematic biases in the model, indicates 
that when the OLBI and the physical model reduce 
random errors, they also tend to reduce the systematic 
error. 

5. The optimal projected baseline lengths for observation 
of P Cyg with a nearly ideal interferometer lie in the 
range 45-90 m. 

6. If several model parameters are to be determined si- 
multaneously, the OLBI observations at a single pro- 
jected baseline cannot provide a substantial error re- 
duction. 
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7. Final remarks 

In this section we discuss to what extent our results are 
relevant for real observations, and how the methods de- 
veloped here could be generalized to more complex cases. 

7.1. Additional sources of measurement errors 

Although general formalism presented in Sect. ^ remains 
valid for any source of errors provided that the linearized 
least-square method is applicable, Eqs. ( pi] ) and ( |l8| ) need 
to be modified if the accuracy of measurements is not shot 
noise limited. 

In the simplest case of uncorrelated errors in spectral 
channels, this modification reduces to multiplication of the 
integrands in Eqs. ( pd] ) and ( p^ by corresponding ratios 
(shot noise error) / (total error) . If this ratios are nearly 
constant across the observed spectral range, then after 
appropriate scaling our results can be directly applied to 
real observations. 

The assumption of uncorrelated errors may be invalid 
for interferometric observations, the observables Wi being 
obtained as a result of much more complicate d proc essing 
of raw observational data (cf. Mourard et al 1994 ) then 
that considered in This point must be further analysed, 
and in future works, it may be necessary to use directly 
Eq. (I) instead of Eq. (|l|). 

In any case, a generalization on arbitrary measure- 
ment errors could be easily incorporated into existing code 
with no considerable increase in the required computer re- 
sources. 



7.2. Observation in several spectral lines 

The theory developed in Sect. || can be applied to the 
case of multiple observed lines without any modifications. 
Formally, it reduces to changing integration domains in 
Eqs. dnl) and (|l|). 

Rapid evaluation of source function, on the contrary, 
would require a detailed study and development of efficient 
methods specific to each spectral line involved, because 
different lines can be emitted in different regions of the 
envelope and originate from different physical processes. 

7.3. Multi-baseline OLBI of spherically symmetric objects 

The method used here can be generalized to visibility mea- 
surements done simultaneously at more than one baseline: 
the only modification required is the replacement of RHSes 
of Eq. ( p^ by the sum of corresponding expressions over 
given set of projected baseline lengths, the directions of 
the baselines being unimportant because of circular sym- 
metry of intensity distribution. This can be done by mod- 
ifying only the co de p erforming the second stage of calcu- 
lations (see Sect. 14) and using the data prepared for the 



single-baseline case during the most time consuming first 
stage of computation. 

7.4. Asymmetric brightness distribution and information 
on fringe phase 

Although time-averaged structure of P Cyg wind can be 
considered spherically symmetric, polarimetric (Taylor et 
al. |199l| ) and interferometric (Vakili et al. 1997) obser- 
vations revealed existence of time-dependent "clumpy" 
structures in the envelope, giving rise to deviation from 
spherical symmetry. 

In this case, measurements of the fringe phase can yield 
important information. This contrasts with the spherically 
symmetric case, where the fringe phase is irrelevant to 
physical modehng of the object. 

On the other hand, a unique determination of model 
parameters using only spatially unresolved spectroscopy 
appears impossible for the complex physical models result- 
ing in two-dimensional brightness distributions, so that 
even the formulation of the comparison problem for spec- 
troscopy and interferometry can be hardly done. The 
choice of the optimal set of baseline lengths for interferom- 
etry still remains the important problem, and the method 
developed here can be used for such an optimization after 
certain modifications. 

First, the fringe phases (or observable linear combina- 
tions thereof) should be included into the set of observ- 
ables defined in Eq. (p7|), and Eq. ( p^ ) should be modified 
accordingly. This would require comprehensive error anal- 
ysis depending on particular method used for extracting 
phase information from measurements. 

Second, errors in parameter determination for obser- 
vations performed with interferometer configurations Bi 
and B2 (here Bi and B2 denote corresponding sets of 
simultaneously used projected baseline lengths) can be 
compared using the ratio 

Ci(0,r,AA,Bi)/Ci(©,T,AA,B2) instead of 
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Appendix A: Computation of DiOj and M{niWi) 

In this Appendix, we compute the non-trivial elements of 
the covariance matrix of errors used in Sect. 3.2. Since we 
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consider here a single spectral channel, the index i will be 
omitted everywhere. Thus, the values N, N, n, W, W, and 
w in this Appendix are identical with the values denoted 
as Ni, iVj, m, Wi, Wi, and Wi respectively in Sect. |[ 

In practice, the error in W depends on too many details 
of experimental techniques, atmospheric conditions and 
data processing methods to be analyzed in general form 
(Mourard et al 1994). The present analysis is restricted to 



a highly idealized situation. Namely, as well as in Sect. 3.1 



the only source of measurement errors taken into account 
is the photon shot noise. We assume that the static fringes 
are detected by counting the photons in nonintersecting 
channels that cover in total the interval of length Z of 
optical path difference (OPD). If z is the OPD, which is 
measured here in units of X/2tt, and F{z) is the photon 
counting rate per unit z, then 



F{z) ^ N{l + Vcos{z + zq))/Z , 



(A.l) 



where zq is unknown fringe phase shift and N is propor- 
tional to the total received flux. 
If we define 



C ^2 [ F{z)coszdz 
Jo 



and 



5 = 2y F{z) sin zdz, 



(A.2) 



(A.3) 



then assuming that Z ^ 1 and neglecting the terms that 
decrease as 1/Z for Z —> oo we obtain that 



N = 
and 



/ F{z)dz 
Jo 



(A.4) 



(A.5) 



An estimate for W — can be cons tructed by replacing 
the values in the RHS of Eq. ( [A.5| ) by t he ir es timates 
that can be easily deduced from Eqs. ( [A.2| )-( A.4). If the 
interval of OPDs covered by measurements is divided in M 
subintervals of length ApZ and ApZ ^ 1 for p = 1, . . . , M, 
then estimate W oiW based upon the number of photons 
received in each of these subintervals is given by 

W ^ {C^ + S^)/N^ , 

where N = J^pLi^pi is the number of photons re- 
ceived in the p-th subinterval, C — 2^^!]^ Np cos Zp, and 



2e; 



M 



Np sin Zp. 



If we assume further that W as a function of Up — Np — 
Np can be linearized in the vicinity of zero and take into 
account that M{np) = and M{npnq) = SpqNp, then after 



some algebra we obtain for the statistical characteristics 
of 



D(w) 




41^ 



i2~W), 



which proves the validity of Eqs. (jl^) and (|l6|) for i = j. 



Appendix B: Simple approximation for the 
escape probability 

In the Sobolev approximation, the escape probability for 
a line photon in the spherically symmetric accelerating 
outflow is given by (See e.g. Mihalas |197S|, Sect. 14.2) 



Jo xi/Qip) 

where 

Qiii) = ^i^dV/dr + (1 - ^i^)iV/r) 



(B.l) 



XI 



(7re^/mc)/ 



V ~ v/vth, V and Wth are the outflow and scattering atoms 
thermal velocities respectively, Aj/^ = j/uth/c, v is the line 
frequency, / is the oscillator strength, and n is the number 
density of atoms on the lower level of the transition. 

Here and in what follows we consider the fixed point 
of the envelope at the distance r from the center and do 
not show explicitly the dependence of (3, V, n, and other 
physical parameters of the envelope on r. Also, in expres- 
sion for xi neglected the stimulated emission, which 
do not introduce a noticeable error when the formation of 
Ha in the P Cyg envelope is considered. 

Our purpose is to find an approximation P to (3 that 
permits fast evaluation. In order that in solving the equa- 
tions of statistical equilibrium the approximation could be 
used instead of exact definition, the former should retain 
principal mathematical properties of the latter refiecting 
underlying physics. Namely, we require that /3(0) — 1 and 
(3{n) should monotonically decrease to zero as n —^ oo. 

The simplest function that obey the above stated re- 
quirements is 



1 + n/us 



(B.2) 



where rias is a parameter. To choose the optimal value of 
rias, we consider the behavior of for large n. When 



n 



oo, the numerator of the integrand in Eq. (B.l 



can 



be replaced by 1 and asymptotic form of (3 can be easily 
computed. 
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Comparison of that asymptotic form with Eq. (B.2) 
shows that if we set 



(71 1 , A dv 2v. 



(B.3) 



then = 1 up to the terms of order l/n when n —>■ oo. 
Here gi and are statistical weights of the lower and 
upper levels of the transition respectively, A is the line 
wavelength, and A^i is the spontaneous emission coeffi- 
cient. 



To study the precision of approximation ( B.2 ) , we per- 
formed numerical computations for a wide range of val- 
ues {dv / dr) / {v / r) and njn^ (it can be easily shown that 
the error is uniquely determined by the values of these 
two parameters). Our results show that the relative error 
reaches its maximum of 0.23 for ^^^^ — that is when 
v(r) oc r, and 7T,2/nas ~ 1-8. 

The approximation found here is applicable if the es- 



cape probability is given by Eq. (B.l), that is if the fre- 
quency change in scattering is described by complete fre- 
quency redistribution. Of course, it is useful only if the 
usual conditions of applicability of the Sobolev approxi- 
mation are satisfied. The expressions (B.2) and (|B.3[) can 



be generalized in two ways. First, the error can be reduced 
if rational approximations of higher order are used instead 
of Eq. (B.2). Second, a similar approximation can be ob- 
tained for a general three-dimensional flow as it will be 
showed in a future work. 
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